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ABSTRACT 

The problem of neutral-current processes (neutrino scattering, pair emission, pair 
absorption, axion emission, etc.) in a nuclear medium can be separated into an 
expression representing the phase space of the weakly interacting probe, and a set of 
dynamic structure functions of the medium. For a non-relativistic medium we reduce 
the description to two structure functions S'^(a;) and Syi^^) of the energy transfer, 
representing the axial-vector and vector interactions. Sy is well determined by the 
single-nucleon approximation while Sa may be dominated by multiply interacting 
nucleons. Unless the shape of Sa{uj) changes dramatically at high densities, scattering 
processes always dominate over pair processes for neutrino transport or the emission 
of right-handed states. Because the emission of right-handed neutrinos and axions is 
controlled by the same medium response functions, a consistent constraint on their 
properties from consideration of supernova cooling should use the same structure 
functions for both neutrino transport and exotic cooling mechanisms. 



I. INTRODUCTION 



The neutrino signal from the supernova (SN) 1987A has confirmed our basic under- 
standing that type II supernovae result from the core collapse of massive stars [1]. 
The observed neutrinos were radiated from the "neutrino-sphere" with a luminosity 
and energy spectrum commensurate with expectations based both on broad theoret- 
ical grounds and detailed numerical models. Further, the time scale for the neutrino 
"light curve" was broadly in agreement with the hypothesis that energy transport 
in the hot, dense core is dominated by neutrino diffusion. Despite these successes 
there is much that needs to be done to achieve both a quantitative and a qualitative 
understanding of SN physics. 

In this paper, we examine a number of issues concerning neutrino transport in 
SN cores. We are especially interested in the neutral-current (NC) processes that 
govern the transport of fi and r neutrinos. These include scattering uX X'v, pair 
emission X — > X'uV, and pair absorption uVX — > X' where X and X' are config- 
urations of one or several particles of the medium. These processes are intimately 
related because the underlying matrix elements are structurally the same. Still, the 
pair processes have usually been neglected in transport calculations [2], presumably 
because in a non-interacting gas of nucleons they vanish due to energy momentum 
conservation. We would like to know if this naive expectation is born out in a strongly 
interacting medium. 

Although this question has always been present, it has been obscured by the 
difficulty of calculating any of the relevant processes reliably. Further, the authors 
that worried about scattering processes and transport in hot SN cores [3,4,5,6], 
could plausibly assume that scattering was the dominant part of neutrino transport; 
whereas authors worried about pair processes were usually concentrating on the cool- 
ing of older, cold, degenerate neutron stars [7]. There is another problem in that 
the relevant calculations are usually done in the context of perturbation theory — 
processes with the fewest number of nucleons are assumed to dominate. The lowest- 
order emission process is usually neutrino pair bremsstrahlung involving two nucleons, 
NN — > NNuu, whereas the lowest order scattering process is scattering from a single 
nucleon uN Nv. Therefore, if pair processes turned out to be important com- 
pared to scattering, the whole perturbative framework of the calculations would seem 
questionable. 

In this light, we find several recent papers on the SN emission of hypothetical 
right-handed (r.h.) neutrinos extremely interesting [8,9,10]. On the basis of per- 
turbative calculations of nucleon-nucleon and pion-nucleon interactions (if there is a 
pion condensate) it appeared that pair processes of the type X — > X'Viv^ would be 
more important than spin-flip scattering u^X — > X'vji in a non-degenerate nuclear 
medium. As these processes are structurally very similar to the ones involving only 
left-handed (l.h.) neutrinos we are led to wonder if pair processes are also important 
for the standard neutrino transport in newborn neutron stars, and under what, if 
any, conditions the perturbative framework can be trusted. Similarly, resolving the 
question of pair processes in transport calculations will affect the interpretation of 
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the emission rates calculated for several species of hypothetical particles, including 
r.h. neutrinos and axions. 

General answers to these questions are difficult because of our lack of knowledge 
about strongly interacting systems. In the end we travel the perturbative road, but 
before doing so we want to make sure that the relation between scattering and pair 
processes is kept clear. To this end, we factorize the problem of neutrino interactions 
into a "neutrino part" , and a "medium part" . This allows us to maintain a consis- 
tent treatment of the medium while discussing the different neutrino processes. The 
medium is then described by a small number of response functions, common to all NC 
processes. In the limit of non-relativistic medium constituents and ignoring the neu- 
trino momentum transfer to the nucleons, the medium response can be reduced to a 
single structure function S{u!) of the energy transfer. 5" is a linear combination of the 
density and spin-density dynamic structure functions at vanishing three-momentum 
transfer. Although we derive S for ordinary l.h. neutrino interactions it also applies 
to spin-flip scattering. 

In a sufficiently dilute medium the dynamic structure functions can be calcu- 
lated perfectly well by perturbative methods, so we may borrow readily from other 
authors' calculations. Our hope is that when extended into the regime of high densi- 
ties the overall shape of S{u!) will not abruptly change. In this case we will argue that 
the relative strength of scattering almost certainly exceeds that of pair processes for 
a medium dominated by non-relativistic nucleons. Thus, if a pion condensate would 
substantially add to X — X'vjiPi it would add even more to X'^R- Then, 

however, it would also strongly affect X'v^ and thus neutrino transport. 

Turn back to the emission of weakly interacting particles. The bremsstrahlung 
emission of r.h. neutrinos or axions are sensitive only to the spin-density fluctuations. 
Even so, the spin-density structure function used for r.h. neutrino or axion emission 
is the same as that for ordinary l.h. neutrino scattering and, therefore, the uncer- 
tainties of, for example, the axion emission rate are not unrelated to that of neutrino 
transport. It follows that a consistent treatment of "exotic" particle emission would 
have to rely on a common structure function S{ijj) for both neutrino transport and 
particle emission, whether or not S{ijj) can be reliably calculated. Modifications of the 
neutrino transport and the emission of novel particles will both affect the observable 
neutrino signal, and thus they should be implemented on the same footing. 

There remains the problem of how to calculate the dynamic structure functions, 
a problem which is arguably more difficult for the present application than for any 
other problem concerning neutrino interactions with stellar material. The medium 
at the core of a proto- neutron star is hot, dense and strongly interacting. There are 
no parameters that one may use as the basis for a perturbation expansion. Further, 
the material is neither degenerate nor non-degenerate. In this context, we suggest 
a phenomenological approach that we believe illustrates some of the features that 
a full treatment of an interacting medium must possess. This approach uses the 
interactions of the medium to regulate soft processes, i.e. those where the energy 
transfer is small, but leaves the hard processes relatively unscathed. Although not 
totally successful, the technique goes in the direction of regulating the amplitude of 
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the response function calculated in a perturbative series, while modifying its shape 
in a controlled way. 

The rest of this paper is organized as follows. In § II we develop the main 
tool for our investigation, the vector and axial structure functions Sv{uj) and 5*^(0;), 
respectively, which describe the medium's response to NCs. In § III wc present useful 
properties of these functions and evaluate them perturbatively. In § IV we examine 
the relative strength of pair vs. scattering processes under a variety of different as- 
sumptions about the high-density behavior of Sa- In § V we consider applications to 
non-standard physics, i.e., the emission of novel particles and the impact of a pion 
condensate. Finally, in § VI, we summarize our discussions and conclusions, with 
some suggestions as to useful strategies for future research. 



II. NEUTRAL-CURRENT NEUTRINO PROCESSES IN A MEDIUM 

a) Collision Integral 

Neutrino transport in a medium is governed by the Boltzmann collision equation, 
L[f] = C[f,f] = {df/dt)coii where / and / arc the neutrino and anti-neutrino oc- 
cupation numbers for l.h. states. (For now we focus on the standard model without 
r.h. neutrinos.) An analogous equation applies to /. The Liouville operator is given 
by L[f] — dtf + X- Vx/ + k- Vk/ while the coUision integral is 



dt 



coll 



+ W^-fe2,fei(l - /ki)(l - /ks) - W^fcl,-fe2/kl7k2 



(2.1) 



where we have written the momentum variables as subscripts. The first term corre- 
sponds to neutrino scatterings into the mode ki from all other modes, the second term 
is scattering out of mode ki into all other modes, the third term is pair production 
with a final-state neutrino ki , and the fourth term is pair absorption of a neutrino of 
momentum ki and an anti-neutrino of any momentum. In this collision integral we 
only include effective NC processes between neutrinos and the nuclear medium, i.e., 
we are not considering charged- current processes and NC processes coupling neutrinos 
to the leptons in the medium. 



^ The temperature T is expected to be of order 50 MeV, whereas the density of 
baryons is characterized by Fermi momenta in the 400—500 MeV range. Particle 
species whose density is characterized by are not nearly as abundant as 
nucleons. The importance of degenerate e~ and z/g for neutrino scattering is 
typically less than a tenth that of the nucleons. Even though it is conceivable 
that there may be situations where they dominate [11], we will always ignore 
purely leptonic processes. 



4 



The expression W^^^j^^ = W{ki,k2) is the rate for a neutrino in state ki to 
scatter into k2 via interaction with the medium. W is defined in the entire k space 
for both of its arguments where negative energies indicate the "crossing" of an initial- 
state neutrino into a final-state anti-neutrino or vice versa. In thermal equilibrium / 
and / must assume Fermi-Dirac distributions at a given temperature T and chemical 
potential ii while the collision integral must vanish, leading to the detailed-balance 
requirement W{ki,k2) e^'^^-'^^y^W{k2,ki). 

It is the pair production and absorption terms under the integral in Eq. (2.1) 
that are usually neglected in the context of SN neutrino transport. The justification 
appears to be that the pair terms vanish for a medium of free nucleons. More specif- 
ically, W{ki, /C2) is identically zero if the energy momentum transfer k — ki — k2 is 
timelike {k^ > 0) and if the dispersion relations of the medium excitations are like 
that of ordinary particles, E"^ — — m? > 0. 



h) The Transition Rate W{ki,k2) 

The low-energy NC Hamiltonian for the neutrino field ^0 and an effective current 
operator i?^ for the medium is 

^int = ^i^'^?^7M(l-75)V'. (2.2) 

If the neutrino interactions are dominated by nucleons, the current of the background 
medium is 

B^^Y. ^a^iCv,i - CA,a5)^i , (2.3) 

i=n,p 

where ipn,p are the interacting quantum fields for protons and neutrons, and Cv,i and 
CA,i are the relevant vector and axial weak charges. 

The current-current structure of ifint allows the transition rate to be written in 
the general form 

W{ki,k2)^^i^S^,m-, (2.4) 

where Gp is the Fermi constant and Ng the number density of baryons. N^^'^ entails 
the neutrino kinematics, and S^i, describes fluctuations in the medium that produce, 
scatter, or absorb neutrinos, i.e., it is the dynamic structure function of the medium 
for NCs. 

Up to a normalization factor, N^^'^ is given by the neutrino part of the squared 
matrix element for the scattering — > [12] 

N''" = TT-^ + k^ki - kik2 9^" - ^e"^'^^A;l«%l . (2.5) 

Then, by the usual crossing relations the expression for the emission of a pair V\^^i'\^^ is 
given by ki —ki, while k2 — > —k2 is for the absorption of i^i^J^i^^. Both operations 
leave N^^^ unchanged because we included a factor {2ijJ\2uj2)~^ in its definition. 
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The dynamic structure factor is then functionally dependent only on the energy 
momentum transfer so that a single tensor Sn^{k) describes all processes in Eq. (2.1). 
It can be written directly in terms of the fluctuations in the weak NCs [13], 

S^'^'ik) = dt e'""^ {B^'it, k)S^(0, -k)) , (2.6) 

where k = ki — k2 = (a;, k). By the expectation value (. . .) of some operator we mean 
a trace over a thermal ensemble of the background medium. 

Equation (2.4) is a first-order perturbative result in the weak Hamiltonian. 
The medium, on the other hand, is strongly interacting, so a perturbative calculation 
of Sfj^iy may or may not be possible. Of course, in principle some insight into its 
properties can be gained by laboratory experiments, e.g. particle emission from the 
hot and dense systems produced in heavy ion collisions. 

However, for practical SN calculations one proceeds with a perturbative ap- 
proach as follows. For a particular term in the expansion one calculates the squared 
matrix element, which can be written in the form \M\'^ = {Gl/8) M^^N^"'2uji2uj2, 
where M^'^ is the square of the medium part of the matrix element. Then one per- 
forms a phase-space integral over all medium participants and sums over all relevant 
processes. 



Nb proSsesH/ 2E,(2^)3 J^j2Ef{2ir ^ ^^^^ 
x{27:)U^ {k + Y.Pi-j:Pf] 

i=l f=l / spins 



(2.7) 



where and Nf are the number of initial and final medium particles, the /p are 
occupation numbers, and the (1 ± /p) are Pauli-blocking or Bose-stimulation factors. 

Unfortunately, for the conditions of a SN core higher-order processes typically 
yield contributions to 5*^1/ which are of the same order or larger than lower-order 
ones (see § III below), so one may question the validity of this method. However, 
because W{ki, /C2) factorizes, this complication can be separated from a discussion of 
the relative importance of scattering vs. pair processes, which is largely an issue of 
neutrino phase space. 

Therefore, instead of relying too much on perturbative calculations of S^'^{k) 
we should focus on its general properties. To this end we assume that the medium 
is homogenous and isotropic, in which case S'^'^ can be constructed only from the 
energy-momentum transfer k and the four- velocity u of the medium [14] , 

S^"" = Ri u^'u''+R2 {u^'u''-g^"')+R3 k^'k''+R4 {k^'u''+u^'k'')+iR5 e^"'°'^Uakp . (2.8) 

The structure functions . . . , i?5 depend on the medium temperature and chemical 
composition, and on the Lorentz scalars that can be constructed from u and k, namely 
A;^ and uk. (The third possibility = 1 is a constant.) Instead of k^ and uk we will 
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use as independent variables the energy and momentum transfer uj and |k|, measured 
in the rest frame of the medium. 

It is convenient to calculate the interaction rate in the rest frame of the medium 
defined hy (1, 0, 0, 0). Multiplying Eqs. (2.5) and (2.8), we find 



GINb 



; 1 + COS e) i?l + (3 - COS ^) i?2 - 2 ( 1 - cos 9) {ui +UJ2)R5 



, (2.9) 



where 6 is the neutrino scattering angle. It is specific to the contraction with N^^, 
relevant for l.h. neutrinos, that only i?i,2,5 contribute. For the spin-fiip processes 
discussed in § V we will find that all i?'s contribute while for axion emission only R2 
survives. 



c) Non-Relativistic and Long- Wavelength Limit 

Within the context of a nucleonic medium there are two closely related limits that 
are often taken — the "non-relativistic" and "the long-wavelength" approximations. 
The justification for both is that the nucleon mass is larger than any other energy 
or momentum scale in the problem. (Even though nucleon Fermi momenta are in 
the 400 — 500 MeV range and the effective nucleon mass may be as low as 600 MeV, 
it is still reasonable to treat the nuclcons non-rclativistically for the purposes of this 
paper.) Then, the currents are expanded as a power scries in l/m^v, after which 
only the leading term is kept.^ The long- wavelength approximation assumes that 
the three-momentum transfer to the nucleons is small compared to typical nucleon 
momenta and thus can be ignored when calculating the available phase-space, thus 
simplifying those calculations considerably. 

As a consequence of the non-relativistic assumption the structure functions 
-^3,4,5 vanish. To see this we first note that the medium current in Eq. (2.2) generally 
is a sum of a vector and an axial- vector piece, = V^+A^^. Then the transformation 
properties under parity of the five terms in Eq. (2.8) arc easily identified to be like 
(V^V) or (A^A^) for terms 1 - 4 and {A^'V'' + VA"") for term 5. In the limit of 
non-relativistic nucleons has only a zero- component, = Cyil'^i^ (a sum over 
nucleon species is implied), while the axial current A'^ has only spatial components, 
A* = Cyi'j/'VV) where a* arc the Pauli matrices and the non-relativistic are Pauli 
spinors. If = and A'^ = all terms involving (V'^V^) or (A'^A^) vanish, i.e., 
terms 3 and 4. Term 5 involves components {A'''V^) which also vanish. Term 1 
only has a 00 component and thus corresponds to {V'^V'^), i.e., for only one species 
of nucleons Ri — CySp with Sp, the usual dynamic structure function for nucleon 
density fluctuations [6] . Term 2 only has spatial components and thus, it corresponds 
to {A'- A''), i.e. B.2 — C\Sfj with S(j the dynamic structure function for spin-density 



^ It may happen that the leading contribution to the current makes no contri- 
bution to an interaction, in which case the next term in the non-relativistic 
expansion should be kept. 
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fluctuations [6]. For a mixed medium of protons and neutrons, tlie interpretation of 
i?i^2 is more complicated because there are isospin and 1 contributions to both [4]. 

In the non-relativistic limit we are thus left with only two structure functions 
Rl,2{k) to consider. Turning to the long- wavelength approximation we assume that 
the three-momentum transferred to the nucleons is negligible. This is justified by 
considering the perturbative series Eq. (2.7). Each contribution to 5*^1/ "knows" about 
the energy-momentum transfer only through the overall energy-momentum conserving 
5-function. If the nucleon mass is very large we may neglect the momenta of the 
neutrinos (and other relativistic participants) in the overall law of three-momentum 
conservation so that 



(2.10) 



It may seem somewhat strange to drop the momentum transfer and keep the energy 
transfer. Specifically, when scattering from individual nucleons in the non-relativistic 
limit, the energy transfer is smaller than the momentum transfer: uo = k-p/m^v < |k|, 
where here p is the nucleon momentum. Although this is true for single-particle 
scattering, and indeed implies that a; = in the long-wavelength approximation, it is 
not true for the higher-order terms involving the interaction of a neutrino with a pair 
(or more) of interacting nucleons. In that case, it is possible to conserve momentum 
exactly in the nucleon sector while still releasing energy to or absorbing energy from 
the neutrinos. Thus, in the long-wavelength limit the perturbative response function 
is a 5-function in cu at leading order, but at higher order the medium can transfer any 
amount of energy subject to the thermal constraints of the medium and the energy 
available in the incident neutrino; but to all orders it is independent of k. 

Thus, in an isotropic medium there remain two structure functions 6*1/(0;) = 
Ri{u!, 0) and Sa{i-o) = i?2('^, 0) which are given in terms of field correlators as 



1 

3Nbj-oo 



Nb J-oo 

SaH = ^ C^te^^' {A{t) ■ A(0)) , 



(2.11) 



where 



V\t) = j d^x 
A\t) = j d^x 



CA,p'4^l{t, yi)(T''ll)p{t, X) + CA,ni^i{t, x)(7Vn(^, x) 



(2.12) 



with the Pauli matrices o"*. 

For an isotropic distribution of neutrinos, after averaging over the scattering or 
emission angles the cos^ terms in Eq. (2.9) average to zero. Therefore, we are left 
with W{ki, k2) — G'yNbS{uj) with a single structure function 



S{u) 



Sy{u) + 35*^1(0;) 



(2.13) 
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In the non-relativistic and long-wavelength limit this function is all we need to study 
neutrino scattering and pair processes. 



III. EVALUATION OF S{uj) 

a) Simple Properties of Sy and Sa 

In general, a full evaluation of S{u!) is not possible. We will present perturbative 
calculations later in this section, but there are hmitations to the range of their reli- 
ability. However, even without detailed knowledge of the structure functions we can 
easily state several useful and simple properties of Sv{uj) and S'^(a;). 

To begin, the detailed-balance requirement for W{ki,k2) translates into 

Sv,a{-^)=^~"^^Sv,aH (3-1) 

so that it is enough to specify these functions for positive energy transfers (energy 
absorbed by the medium). 

Next, we consider the overall normalization. By integrating Eq. (2.11) over 
doj /27r we extract the total strength. On the r.h.s. / e*'^^ du; /27r — 5{t), allowing us to 
trivially perform the time integration, and find 

J-OO ZTT IMb (o n\ 

r+°o du 1 ^ ' 

Of course, we could have chosen any time, not just t — 0, as a, reference point. 

If there are no spin or isospin correlations between different nucleons of the 
medium we may interpret A as the single-nucleon spin operator (apart from overall 
factors). In this case Eq. (3.2) yields the "sum rule" 

/ Tr^vil^) = XpCyp + XnCy^n = '^V , 

(3 3) 

where Xp and Xn are the number fractions of protons and neutrons, respectively.^ 



^ If the medium is degenerate, but approximated by non-interacting particles, one 
must account for final state blocking effects and replace the number fractions 
by Xi — > N^^ f d^p 2/j(l — /j) / (27r)^. This result, however, cannot be extended 
trivially to interacting nucleons. 
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For non-interacting nucleons the time evolution of the operators is trivial, e.g. 
A{t) = A(0), and so we find 

Sv,AH^Sv,A'2n5{u;). (3.4) 

This result yields directly the usual neutrino scattering rate in a medium of non- 
relativistic, non-degenerate, qiiasi-free nucleons: F = G^Nbuji{Sv + 3<Sa)/47r, where 
oui is the energy of the incident neutrino. 

It is expected that the normalization and shape of the response functions will be 
altered in an interacting medium, however this primarily affects the axial response Sa- 
For the vector current Eq. (3.4) should be a good approximation even in the presence 
of nucleon-nucleon interactions because the quantity V*^ (the zeroth component of 
the four velocity) does not fluctuate; i. e. it is identically 1 in the non-rclativistic limit 
where the wavelength of the momentum exchange is large compared to the smearing 
of the nuclear charge density due to interactions in the medium. The axial charge, on 
the other hand, is a spatial 3-vector quantity. Although a nuclcon's location is not 
smeared out, nucleon interactions in the medium cause the direction of the nuclear 
spin to fluctuate, altering the effective charge seen by a leptonic probe. We expect 
that the axial response is weakened for soft energy transfers, cu < Fint, where by Fi^t 
we mean some typical rate for the nuclear spin to fluctuate due to collisions.^ 

The counterpoint to these arguments is that there is no vector response for 
non-zero energy transfers, but in general the axial response gets spread out to higher 
ui with increasing interactions. Specifically, explicit calculations (see § lll.b below) 
show that only the axial- vector current contributes to bremsstrahlung in the non- 
relativistic limit. 

A related issue concerns the appropriate values to use for the axial and vector 
charges Ca and Cy- In a dilute medium, one should use the vacuum values, but 
in a nuclear medium the axial weak coupling constants may be altered. The vector 
charges remain unchanged. Specific values are discussed in Appendix A. 

Thus, within the non-relativistic and long-wavelength approximation the vector 
current is unaffected by multiple collisions but docs not contribute to bremsstrahlung, 
and the reverse applies to the axial- vector current. It remains to determine the axial- 
vector structure function 5'yi(a;) in a fully interacting medium. 

Finally, it should be clear that random spin fiuctuations of the sort which 
happen in a thermal medium may lessen the total scattering interaction rate, but 
we would be surprised if any significant enhancement could occur through thermal 
interactions. Thus, the normalization of S in Eq. (3.3) should constitute a rough 



t Isospin fiuctuations do not lead to an equivalent damping of the vector response, 
due to isospin conservation in the strong interactions. In order to have the 
equivalent effect a proton must turn into a neutron, but this can only occur if 
a neutron turns into a proton in the same interaction. In the long wavelength 
approximation, the net weak vector charge probed by the neutrinos will remain 
unchanged. 
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upper bound. This does not allow, however, for collective effects such as spin-waves, 
or the influence of a pion condensate, which could in principle enhance the interaction 
rates of neutrinos with the medium. 



b) Perturbative Result for S{u!) 

We proceed with a perturbative calculation of Sa{oo). After the 5-function con- 
tribution from free nucleons, the next lowest-order process that can contribute is 
nucleon-nucleon scattering. We consider non-degenerate nucleons and model the in- 
teraction by one-pion exchange. Neglecting ttIt^ in the pion propagator we extract 
from Friman and Maxwell's matrix element [7] and from Turner's emission rate for 
the bremsstrahlung process NN — > NNVv [8] the result 

SNN{uJ) = C\^s{uJ/T). (3.5) 



where the NC axial coupling Ca is, for simplicity, taken to be the same for neutrons 
and protons, and 

r^^/(X„)48.V2 4iV^_ ^3_g^ 

Further, f{Xn) = 1+4X^(1 - X„), and a^r = {f2mN/m-,^)'^ /At: 15, where / fa 1 is 
the pion-nucleon coupling. The function 

s{x) = dy (]f + xyj e"^ !^ M + - a;j (3.7) 

with X — UJ /T is the remnant of the nuclear matrix element after doing the phase 
space integration for the nucleons. The analytic approximation has a maximum error 
of less than 2.2 %. 

This result cannot be the complete answer, however, as it diverges at = 0. The 
perturbative calculation is based on the assumption that a given bremsstrahlung event 
can be viewed as a single collision where in-states travel unperturbed from the infinite 
past and are scattered into out-states which remain unperturbed into the infinite 
future. In a medium these assumptions are violated as each nucleon participant in a 
given reaction has emerged from a previous interaction and will interact again in the 
future. The nucleon field can freely evolve only for an approximate duration r = 
where Fint is a typical nucleon-nucleon scattering rate. Hence, we cannot expect such 
a calculation to yield meaningful results for energy scales = Fi^t. This 

estimate reveals the approximate range of uj where a naive perturbative approach can 
be trusted. 

In a non-degenerate medium in thermal equilibrium, the scale of a typical energy 
transfer in any collision or radiation event is given by the temperature T. Hence, as 
long as T/Fint S> 1 one may expect that a naive application of pertubative methods 
is justified and we call such a medium "dilute" for the purposes of our discussion. 
Conversely, the criteria T/Fint ^ 1 quantifies our notion of a "dense" medium (the 
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limit of large collision rates) where the naive use of perturbative results is problematic 
because it is not justified in any obvious way. The hot material of a young SN core 
appears to be a medium which is approximately non-degenerate and yet "dense" in 
the spirit of this definition. 

Energy conservation in the collisions between neutrinos and nucleons will be 
satisfied only with a precision AE pa Tint- This observation suggests that the nucleons 
may be described as resonances with an approximate width Tint, rather than as states 
with a fixed energy [15]. In a previous paper we have argued [16] that in this case 
one should replace the soft behavior of S'^(cj) by a Lorentzian l/{uj'^ + Tf^^/A). 
If we use this form for the soft part of Sa{i^), and identify Tint = ^^A, we may write 

S'NNi^) = C\ s{u/T) . (3.8) 

This form has two attractive features in the limit of a dilute, but interacting, medium. 
For small values of it approaches the appropriate 5-function for Sa while for 
cu ^ Fa the correct bremsstrahlung emissivities obtain. 

c) Large Collision Rates 

We have argued that in the limit of small collision rates, defined by the condition 
F^ ^ T, the low-o; behavior can be reasonably approximated by a Lorentzian shape. 
In a supernova core, however, wc are confronted with the opposite case of large 
collision rates. Numerically, Eq. (3.6) becomes 

££.415i^f^^V'V55^y'7/<^l, (3.9) 
T po \ ruN J \ T J \ 1.84 J ' ^ ^ 

where Po = 3xl0^^g/cm^ is the nuclear density, and wc have scaled the factor /(X„) 
to its value at X„ = 0.3. The numerical value of F^/T is approximately a factor 
10 larger than the value we estimated for Fint/T" in a previous paper [17]. That 
estimate was based on low energy p-p and p-D scattering data, which may or may 
not be appropriate as a model for the source of spin fluctuations in a medium. It is 
the enormous magnitude of F^/T that makes the interpretation of the perturbative 
results for Sa very problematic. 

If wc used the form Eq. (3.8) then both neutrino scattering and pair processes 
would be substantially suppressed. This will be shown explicitly in § IV. For now, it 
suffices to note that in a thermal environment, typical values of uj are a few T, and if 
Ta^ T then the denominator in Eq. (3.8) becomes large for reasonable uj. Although 
there may be scattering strength at large cu that strength will not be sampled by 
thermal energy transfers. 

It is plausible, however, that in a dense medium SaI^^) is narrower than indi- 
cated by a naive extrapolation of the dilute-medium result. The width of the spin-spin 
correlation function for a single nucleon is determined essentially by scattering off the 
spatial spin fluctuations in the medium. But if those spins are all fluctuating rapidly 
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in time, then their spatial fluctuations will be damped by the same mechanism. To 
get a consistent picture, one must self-consistently include the effects of interactions 
in an evaluation of the interaction rates. 

Even though it appears impossible to extract reliable results for S from a simple- 
minded perturbative calculation, we believe that the overall shape of the structure 
function will be a broad distribution, and that the high-a; wings will be represented 
by the bremsstrahlung result. However, it is in no way evident what one should 
use for the soft part of the axial response function. The Lorentzian shape espoused 
in Eq. (3.8) is possible, but it is also possible that some radically different shape is 
correct. Indeed, in the next section we will show that extreme changes to the shape 
arc required if pair processes are to dominate over scattering processes. 



IV. SCATTERING VS. PAIR PROCESSES 



There are many ways that one could choose to compare the strength of the scattering 
and pair processes. For example, we could concentrate on the contribution to various 
transport coefficients or calculate the Rosseland mean opacity. To illustrate our 
arguments we will use the rate at which particles are absorbed by either the pair 
or scattering terms, averaged over the neutrino phase space. 

These quantities have the advantage of being simple to calculate and interpret, with- 
out being too specific to the details in some corner of phase space. We have normal- 
ized by the number density of neutrinos Ni, — J /(a;i)d^ki/(27r)^. The occupation 
numbers will be taken to be Fermi-Dirac distributions without a neutrino chemical 
potential. 

In § II we argued that for the conditions relevant for a SN core W{ki,k2) 
depends only on one function S of the energy transfer u according to Eq. (2.13). 
Therefore, 



scat 



GINb 



P rOO rOO 

^J^ duji du2ululf{ui)[l-f{u2)\S{u), (4.2) 



or transforming to an integral over uj — UJ2 — oj\ — —uj 



scat 



dc^ dcuiLul + ujif fiuji) [1 - f{u+uji)] S{-ij) . (4.3) 
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By means of the detailed-balance requirement S{—uj) = e ^^"^ S{uj), and relabeling 
the integration variable uj ^ ui, we may write this as an integral over positive u) only, 



scat 



fOO 

I du;F^^t(u;)S(u;). (4.4) 



We have arranged things so that we will need 5'(c<;) only for positive energy transfers 
(energy absorbed by the medium). The detailed-balance Boltzmann factor e'^/"^, rel- 
evant for negative energy transfers (energy given to the neutrinos) has been included 
in the definition of the phase-space factor F{u!) which is specific to a given process 
and chosen type of thermal average, e.g. by number, by energy, etc. 

Fgcat is then expressed as a sum of two terms because we split the cu integration 
to ensure that uj is always positive. 



i^scat(t^) = ^^4^ 1 duiufiui-u)^ f{ui) [1 - f {ui-u}] 



+ 



(4.5) 



Note that in the first term uj = uji — UJ2, while in the second uj = 002 — wi- Also, 
there is a detailed-balance factor e~^^'^ in the second integral relevant for the case 
when the medium gives energy to the neutrinos rather than absorbing it. Finally, in 
the case S{u;) ~ 6{u;) only one of the terms in Eq. (4.5) should be included, to avoid 
double counting. 

Similar manipulations lead to Fpair = dcu Fpg^iT.{u!)S{u!) with 

Fpairiuj) = ^ dwi UjI{uJ - UJif f{ui) J{u-ui) . (4.6) 

Because we consider thermal neutrino distributions, Fpair equally applies to the phase- 
space averaged pair emittance. 

As a further simplification, for non-degenerate neutrinos we may approximate 
the Fermi-Dirac by Boltzmann distributions, and ignore the Pauli-blo eking factor in 
Eq. (4.1). In 

-Fgcat ^-iid -^pair this means that f{ui) e-'^i/^ and f{uj2) ^ e-^^^/r 
while 1 — f{uJ2) — 1- Then we find analytically. 



Fscat(x) ^ ^^^^ (24 + 12a; + 2x^) e"- , 
GINbT^x^ _^ 

^pairW- 27r2 60 ^ ' 



(4.7) 



where x = u/T. These approximations are shown in Fig. 1; they deviate from 
Eqs. (4.5) and (4.6) by less than 10%, except for Fpair at values of x less than 5, 
where Eq. (4.7) overestimates Fpair by up to a factor of 3. Since Fpair <^ Fgcat for 
these x, we can use Eq. (4.7) without reservation. 
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From Fig. 1 it is clear that scattering is more sensitive to the structure function 
near a; = while pair absorption is more sensitive to S at higher cj. Put another way, 
for pair processes to dominate, S{u!) would need to have much more power at large 
uj than near uj = 0. For example, if S{uj) would vary like a power law a;", Fpair would 
exceed Fgcat for n > 3, but would be subdominant for smaller n. 

For a non-interacting medium one can use the 5-function result Eq. (3.4). This 
limit allows for no pair processes, and scattering dominates. The scattering rates due 
to this choice of S{uj) are shown as solid horizontal lines in Fig. 2. Including inter- 
actions broadens the axial response. For a dilute interacting medium, the Lorcntzian 
shape, Eq. (3.8) seems a reasonable approximation. For large uj, the axial response 
then behaves as u>~'^, a result which is much softer than the n > 3 power law re- 
quired for pair processes to dominate. Again, there is no scope for pair processes to 
dominate. 

In the limit of large collision rates there is more room for speculation. In order 
to avoid the conclusion that Fgcat > Tpair we would have to assume that 5* (a;) is not 
monotonically decreasing, but rather is "hollowed out" near a; = 0. A radical example 
would be to take Eq. (3.5) for S{uj) above some cutoff uj^nm, and to set S{u) = for 
UJ < cUmin-^ Under this assumption, the rates Fgcat and Fpair ai'e shown as the dotted 
curves in Fig. 2. As can be seen, Snn becomes more important for pair processes than 
for scattering for uj^am/T ^ 5. However, even if such a radical transformation of the 
axial response occurs, there will still be the vector contribution, S{uj) — 27rSvS{uj), 
to scattering. From Fig. 2, for p = pQ and uJram/T ^ 6, vector current scattering 
dominates over pair processes. Thus, even with the radical assumption of a total 
cutoff of the low-cu response there is almost no range where pair processes dominate. 

Not only must the axial response be hollowed out, but in order for the pair 
processes to dominate, the axial strength at high uj must remain strong. At the 
same time, we note that for uj^ij^/T ^ 4, the scattering response is enhanced over 
the single nucleon result, a consequence of the a;~^ behavior of S^n. If we use the 
Lorentzian form, the uj~^ behavior is regulated, but the strength of the axial scattering 
is severely suppressed. For example, using Eq. (3.8), for uj > ujram instead of the naive 
bremsstrahlung result, Eq. (3.5), yields the two dashed curves in Fig. 2. Again, pair 
processes exceed scattering for a cutoff uJmm/T ^ 5, however now both processes are 
two orders of magnitude smaller than the vector contribution. 

The results shown in Fig. 2 are for p — pQ, and F^ 415. As the density 
is raised even higher, the single particle contributions to the rates increase linearly 
with density, while the naive bremsstrahlung rates increase as p^. So it would appear 
that eventually pair processes could overtake the contribution from Sy- However, if 
we look at the contribution from S'j^j^ instead of Snn we find that the contribution 
to F is independent of p, and the comparison favors scattering even more strongly. 



^ For the sake of this extreme example we ignore the normalization of Eq. (3.3). 
If that relation is indeed an upper bound then the implied cutoff is a;min ^ 10^. 
Such a large value suggests that some other form of regulation must apply, e.g. 
as in the Lorentzian form S'j^j^. 
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Given the wide range of possibilities for S, it is not clear whether or not increasing 
the density can eventually allow pair processes to win out. 

We do not think that an extreme behavior such as a sharp cutoff to S{u;) is 
realistic, so we view these exercises as examples of how radically different the high- 
density behavior of S would have to be in order that pair processes might dominate. 

We conclude that, even though we cannot calculate the exact form of the axial 
response of the medium, a combination of phase space arguments and a recognition 
of the stability of the vector response, seems to preclude pair processes dominating 
over scattering processes in the evaluation of neutrino transport properties. In order 
to avoid this conclusion the axial response function for a strongly interacting medium 
would have to be both suppressed at low and enhanced at high energy transfers in 
comparison to the bremsstrahlung result, Snn. 

Another effect of increasing density is degeneracy of the nucleons. As the nucle- 
ons become degenerate the amount of nucleon phase space available for small energy 
transfers decreases, effectively creating a harder response function. From Friman and 
Maxwell's [7] bremsstrahlung rate we can easily extract Snn- The result includes 
various prefactors, but for our purposes the interesting part is the dependence on the 
energy transfer, 

«deg(^) = ^2^:^ ^ , (4.8) 

where x = uj/T and 7 = Fint/T is our phenomenological cutoff, again necessary to 
regulate an divergence at small u. Now, however, the large x behavior is much 
"stiffer" , Sdeg oc x. This is still not hard enough to provide for the dominance of pair 
processes, but it is getting closer. In the extreme case where Fint ^ T, the effective 
power law of the response function is n = 3. This last value is enough to make pair 
processes competitive; however, we don't think this extreme case really applies to the 
SN core since the temperatures are high enough that the nucleons are only mildly 
degenerate. 

In summary, we repeat that it is in no way evident how one ought to extrapolate 
the dilute-medium approximation for SaI^^) into the high-density regime. However, 
unless the overall shape of S is radically different from its limiting behavior, the pair- 
processes will always be less important than scattering. Further, it appears that both 
the pair-rates and the axial- vector scattering rates will be suppressed relative to their 
"naive" values. 

Finally, one may wonder if these conclusions are not in conflict with the rea- 
sonable agreement between the SN 1987A neutrino observations and the expected 
signal duration. However, as the vector-current appears to remain unsuppressed, the 
scattering rates would be reduced only to about 1/5 of their naive values, even if the 
axial-vector contributions were entirely suppressed. It is an interesting question if, 
given the freedom to adjust other parameters, such a reduction in neutrino opacities 
can be excluded on the basis of the detected events. 
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V. NON-STANDARD PHYSICS 



We now turn to exploring how our conclusions affect various arguments related to 
certain non-standard aspects of SN physics. First, we consider the emission of hy- 
pothetical particles weakly coupled to the NC of the medium. This process could 
provide an anomalous cooling mechanism for the core, which in turn could result in 
a diminished neutrino signal of SN 1987A, thus allowing one to constrain the prop- 
erties of the particles in question [18]. In fact, trying to understand the consistency 
of such arguments was the original motivation for this work. Second, we look at the 
possible impact of a pion condensate for both neutrino transport and hypothetical 
particle emission. All these issues are connected through the same medium structure 
functions that control neutrino transport. It is therefore possible, and necessary, to 
treat the novel effects and neutrino transport in a consistent fashion. 



As a first case we consider the possibility that neutrinos have a small Dirac mass rui,, 
which allows for the production of "wrong-helicity" states in the deep interior of a 
SN core [17]. If rrii, is small enough these states can escape freely and thus provide 
an anomalous sink for the heat in the core. If rrii/ is not too small then that heat 
sink would have had observable effects on the SN 1987A neutrino signature, thus con- 
straining nil,. Several estimates of these bounds have been made, often with different 
treatments of the processes which contribute to the emission [8,9,10,11,12,17,19]. 
Specifically, under certain assumptions, the strength of the pair emission processes 
X X'Vj^uji has been shown to exceed the rate for spin-fiip scattering from single 
nucleons vj^N — > N'l/fi [8, 9, 10]. In analogy to our remarks about the importance of 
pair processes in the transport of ordinary l.h. neutrinos; we argue that unless the 
full response is very "hard" spin-fiip scattering still dominates over pair emission. 

We take rUi, to be very small compared with typical neutrino energies. There- 
fore, it is not necessary to distinguish carefully between hehcity and chirahty, and 
so we shall always refer to wrong-helicity neutrinos or anti-neutrinos as right-handed 
and to the corrcct-helicity states as left-handed. Moreover, if rrii, is sufficiently small, 
r.h. states will not be trapped so that their occupation numbers are sub-thermal. 
The energy-loss rate per unit volume due to the emission of r.h. neutrinos is then 

Qvr — QscsX + <5pair with 



where VFfc^^fc^ is the transition probability for the scattering process ^'l(^l) -\- X ^ 
X' v^fyhji) while W-kL,kR refers to X ^ X' + VLiki) + i^RikR). For the energy 



a) Processes With "Flipped" Neutrinos 




(5.1) 
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loss due to uji one obtains an analogous expression; for non-degenerate we have 

Qvr = Qvr- 

In order to determine Wk^^k^ we note that it can be written in the form Eq. (2.4) 
where the structure function Sfj_y remains unchanged. For interactions of neutrinos 
with specified helicities, Gaemers, Gandhi, and Lattimcr [12] showed that the expres- 
sion for A^^^ remains of the form Eq. (2.5) if one substitutes ki — > ^{ki ±mjySj), i = 1 
or 2, where the plus-sign refers to v and the minus-sign to F, and s is the covariant 
spin-vector. For relativistic neutrinos we may consider a non-covariant lowest-order 
expansion in terms of rrii,. Then ki remains unchanged for l.h. states while 



ki ^ (oji, ki) ki = {mv/2ixn) (cuj, -kj) (5.2) 

for r.h. ones. After this substitution has been performed, all further effects of a non- 
zero mv are of higher order so that one may neglect r?v everywhere except in the 
global spin-flip factor.''" Therefore, the neutrino tensor N^^ is found from Eq. (2.5) 
by multiplication with (mi//2cjjj)^, inserting ki = ki, and k2 = {ujr, — k^^). 

Next, wc explicitly contract N'^'^ with S^i/ in the form Eq. (2.8) in order to 
derive the spin-flip equivalent of Eq. (2.9), 



— , , G^Nb ( mu\ 



1 -cos^) Ri + (3 + cos^) i?2 

\^UJRy L 

+ 4a;|(l - cos 9) R3 - 4.u;r(1 - cos 6) R4 + 2{ujr -ujl){'^ + cos 9) R5 



(5.3) 

Following Gaemers, Gandhi, and Lattimer [12] we emphasize that this expression 
differs in more than the factor (rni,/2ujR)^ from the non-flip rate Eq. (2.9). This 
difference is related to the changed angular momentum budget of reactions with spin- 
flipped neutrinos. 

In the non-relativistic limit, however, the contributions of -R3,4,5 can be ne- 
glected as discussed in § II. c above. Moreover, in an isotropic medium the cos^ terms 
average to zero. Then, the spin-flip factOT is the only modiflcation necessary to deal 
with r.h. neutrinos, and we may write W{ki:^,kR) — GYNB{mi,/2ujR)'^ S{uj), which 
involves the same structure function as for the non-flip case. Then, in analogy to 
Eq. (4.4), the energy-loss rates can be written as 

Qi^ duFi{u)S{u), (5.4) 



^ The dispersion relation of neutrinos in a SN differs markedly from the vacuum 
form; in the core the "effective m,^" is several lOkeV. However, m,^ in equa- 
tion (5.2) is the vacuum mass which couples l.h. and r.h. states and thus leads to 
spin flip while the medium-induced "mass" only affects the dispersion relation 
of l.h. states. This view is supported by a detailed study of Pantaleone [20] . Of 
course, for non-relativistic neutrinos the situation is more complicated because 
an approximate identiflcation of helicity with chirality is not possible. 
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where, i stands for "pair" or "scat" . If for non-dcgcncrate neutrinos we use a Boltz- 
mann distribution and neglect Pauli-blocking in Eq. (5.1), the phase-space functions 
are 

Fscat(x) = ^^i^^^(x2 + 6x + 12) e-, 

AT rpA 4 (^■^) 

where amp = 0^7X1^/ At: and x — uj/T. 

A comparison of Eqs. (5.5) and (4.7) shows that the r.h. neutrino emission bears 
a resemblance to l.h. neutrino absorption: in order for pair processes to dominate, 
the response function must be very hard. In § III.c we characterized the "hardness" 
of 5 by a power law a;"^, and noted that the "critical exponent" where pair processes 
were as important as scattering was n — Similarly, for r.h. neutrino emission, the 
critical exponent is n fa 2.57. 

Since the critical spectrum for spin-flip pair emission is slightly softer than for 
l.h. pair absorption, it is slightly easier for pair-proccsscs to dominate, than in the 
case of l.h. neutrino transport. In Fig. 3 we show the r.h. emission results analagous 
to the l.h. scattering results shown in Fig. 2. The emission rates are scaled to Qq — 
Cflip-^B7'^/(47r^). The results are similar to those presented for the transport of 
l.h. neutrinos. The cross-over to pair processes dominating the axial emission occurs 
at uJram./T ~ 4, but a combination of axial plus vector scattering almost always 
dominates even when no attempt (other than the cutoff at a;iiiin) is made to regulate 
the soft uj response. When the Lorentzian form for Sa is used, vector scattering 
dominates at nuclear density. Therefore, by any indication that we can obtain from 
extrapolating perturbative calculations into the high-density regime, the scattering 
processes are dominant for spin-flip emission of r.h. states, as well as for the transport 
of l.h. neutrinos. 

Even though we can be reasonably sure that scattering processes dominate the 
emission rate of r.h. neutrinos there is still uncertainty about the total emission rate 
due to the uncertainty about the axial contribution. If the axial contribution is much 
suppressed then Qsc&ty should give a good lower bound to the emission rate, whereas 
if the axial part is enhanced the emission rate may be much larger than Qscaty- 
We feel that the latter possibility is remote but cannot be ruled out on the basis of 
perturbative calculations alone. 

Even though perturbative calculations do not seem to allow for a reliable cal- 
culation of the absolute magnitude of neutrino scattering rates, the NC transport 
of ordinary l.h. neutrinos and the spin-flip emission of r.h. Dirac-mass neutrinos are 
linked to each other through a single response function. In fact, the functions i^scat(^) 
and Ff^cat (.T") in Eqs. (4.7) and (5.5) even have identical shapes.^ Because modifications 



t If one of Un or Uj- had a large enough mass so that spin-flip processes would be 
important at all, and if it had a not-too-small mixing angle with i/g, a degenerate 
sea of this flavor would be populated [8, 21]. For such a degenerate flavor the 
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of the transport coefficients and of the spin-flip rates would both affect the observ- 
able neutrino light curve, a consistent attempt to constrain Dirac masses would then 
depend on letting the NC transport coefficients and the spin-flip emission rate "float" 
together. 

b) Axions 

The emission of axions is another interesting possibility for an anomalous energy sink 
in the inner core of a SN. Again, the expected impact on the SN 1987A neutrino 
signature was used to constrain the coupling strength and then indirectly the mass 
of these pseudo-scalar bosons [22] . They would couple to the medium according to 

H,^t^^B^d^<f>, (5.6) 

fa 

where (f) is the axion field and fa is the axion decay constant which has dimensions 
of energy. The medium current is presumably dominated by protons and neutrons 
according to 

where the Cn,p are model-dependent dimcnsionlcss coupling constants of order unity. 
The structure of this interaction is of the current- current type so that we may apply 
an analysis similar to that for neutrino processes, except that we replace A^^*^ by its 
axionic equivalent ^'^'^ = k'^k^/2uja, where ka is the four momentum of the axion and 
Ua its energy. Note that by our definition of the momentum transfer k = —ka- 
li axions are weakly enough coupled they are not trapped in a SN core, so that 
their occupation numbers are sub-thermal. In this case the energy-loss rate per unit 
volume is 

The response function S^^'^ is as defined in Eq. (2.6) with the above current B^. In 
the non-relativistic and long-wavelength limit we have 

= 7%2 rduu^e-'^/^Saiu;) , (5.9) 

where Sa = R2,a- Apart from normalization the phase-space function for axion emis- 
sion, cj^e^'^/"^, is identical to that for spin-fiip pair emission. Note that Ri^a vanishes 
because the current is purely axial. Since the axion couplings to the nucleons 
are generally different from the axial current couplings of the nucleons in the weak 

phase-space functions F^{lj) and F^{uj) for ordinary and spin-fiip scattering 
would be different because the former would involve neutrino Pauli-blo eking 
factors for the final-state v^. The two processes would then involve different 
integrals over the common function S{uj). 
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interactions, the isospin structure of i?2,a will differ in detail from that of R2 for 
neutrinos. However, the general considerations behind a calculation of R2 are en- 
tirely analagous in the two cases. Specifically, perturbative calculations done using 
the one-pion exchange approximation show different functional dependences on the 
proton and neutron fractions, i.e. the analogs to f{Xn) in Eq. (4.8), but the same sort 
of soft divergence in the energy transfer a;~^. Therefore, we may use the neutrino 
response function as a template to discuss the properties of Sa{u!), without 

worrying too much about the numerical accuracy of the prefactor. 

From our previous discussions several points should be evident. First, naive 
axion brcmsstrahlung calculations rely on an Sa similar in form to the iSjVAf given in 
Eq. (3.5). Although there is no divergence in the emission rate, that is only because 
the phase space factor has sufficient powers of cu to regulate the divergence. However, 
since that divergence does show up in neutrino transport, it must be regulated in some 
fashion; e.g., the modification to the Lorentzian form in Eq. (3.8), or the sharp cutoff 
at (^min used in Figs. 2 and 3. It seems likely that whatever mechanism is chosen to 
regulate the divergence at small cu, that regulation will result in a suppresion of the 
axion emission rates compared to the naive brcmsstrahlung calculations, perhaps by 
many orders of magnitude. For example, the effect of the Lorentzian regulation is 
between 3 and 4 orders of magnitude at nuclear density, as can be seen by comparing 
the curves labeled Qpair,iViV and Qp^ir NN -^^S- 3. Conversely, if one were to insist 
that the naive calculation of the axion emission rates were correct, one would have 
to conclude that the neutrino scattering rates are much enhanced, in violation of 
our expectation that thermal spin fluctuations can only decrease the net interaction 
rate. To probe the effects axion emission would have on the neutrino cooling curves 
of supernova cores, one must treat axion emission and neutrino transport using the 
same assumptions about the axial response functions of the medium. 



c) Pion- Induced Processes 

It is also interesting to consider non-standard effects in the nuclear medium that could 
grossly affect SN cooling. In particular, the conditions are close to forming a pion 
condensate, and in some equations of state the number of 7r~ can be much in excess 
of a thermal distribution [23]. A pion condensate can have an important impact 
on the late-time cooling of neutron stars [24], and it was recently speculated that 
under certain assumptions it could have a substantial effect on neutrino processes in 
a newborn neutron star also [8. 10]. 

We calculate the one-pion contribution to the structure function on the basis 
of n~p in the hadronic initial state, and n in the final state, while we neglect 7r° and 
tt'^ as we are mostly interested in the case of a 7r~ condensate. For the pion-nucleon 
interaction we use the pseudoscalar form 

Hnnp = V2 V'„757r"'0p ■ (5-10) 
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We set rUn = rup = uin, and keep = (2/mAr/m7r)^/47r constant. Moreover, we 
calculate in the non-relativistic and long- wavelength limit, leading to 5^(kp — k^) 
for momentum conservation. Thus, — Ep and we are left with 5{u!t^ + 00) for 
energy conservation. With these simplifications we find for the one 7r~ contribution 
to Eq. (2.7)t 

^ "'^ J (2^^^^P^ - ^-^P)] J (2^ ^ WTkwW^ ^ 



^2 



(5.11) 

where P'^'^ = g^^ — u^u'^ . The isospin couplings are defined by Ca,o = 2i^Ap~^An) 
and Ca,i = \{Ca,p + Ca^u), and similarly for the vector couplings. The isospin 1 
vector coupling does not contribute in the non-relativistic approximation. Also, we 
have dropped a term proportional to Ca,iCv,q which is formally of the same order of 
magnitude but vanishes when averaged over neutrino directions since it behaves as 
(ki X k2) • k^. 

The term [{p + ■ /c^r]^ in the denominator of Eq. (5.11) comes from the 

propagator of the intermediate nucleon. It differs slightly between the two amplitudes 
which contribute to Eq. (5.11), but we dropped the neutrino 4-momentum from the 
nucleon propagator so that the two diagrams would be more similar. In the non- 
relativistic approximation and for thermal pions, it is equal to (mjvcuTr)^; then the 
integral over the nucleon phase space gives 1. However, for a pion condensate we will 
have to be more careful as this denominator can diverge. 

Note that the tensor structure of Si^^ reduces to the general form of Eq. (2.8) if 
the pions are distributed isotropically. However, if the pion condensate has a non-zero 
momentum, the covariant Lorentz structure will involve the condensate momentum 
and so, the number of distinct response functions increases. 

As a first specific case we discuss thermal pions. In order to determine the 
contribution St^{u}) to our general response function we contract S!^^ with N^^ , and 
integrate over neutrino angles. After this, any angular dependence on kyr has disap- 
peared so that the pion phase space integration can now be done just over uj^. For 



t In deriving Eq. (5.11) we neglected a diagram where the Z boson interacts with 
the 7r~ before it is absorbed by the nucleon. We also neglected A degrees of 
freedom and made no attempt to modify the tt properties for medium effects. 
These issues are discussed in Migdal et al. [24] in the context of a degenerate 
nuclear medium. Although there is no real justification for neglecting the same 
effects here, it does reduce the complexity of the result and allows a simple 
comparison to other work [8]. 
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non- degenerate nucleons the result is, 




X J 

^Xn[l + U-{\^ 

(5.12) 

If one takes = and Ca,i = 0, this agrees with Turner's result [8]. Note that 
the term with cu > ttIt^ corresponds to n in the initial and pTr~ in the final state, 
i.e., to the creation of a pion. In thermal and chemical equilibrium one must satisfy 
A^TT = /^n — so that Xpf^-{\u}\) — e~l'^l/^Xfj[l + /^-(|a;|)], and detailed balance is 
satisfied. 

Since S'^r oc 0(|ci;| — rriT^), it is an explicit example of a "hard" S with all power 
above some threshold, so it is interesting to compare the scattering and pair process 
strength as a function of my^. Evaluating Ti = Fi{uj)STY{uj) duo, with Fi given in 
Eq. (4.7) for i = scat or pair, and using Eq. (5.12) we find that rpair(iC7,-) = rscat(3:^7r) 
for = m-j^/T = 3.9. At this point the scattering processes are suppressed by about 
a factor of 6 from their maximum if we had ignored the pion mass, i.e., rscat(3.9) = 
l/6rscat(0); but the pair processes which are sensitive to higher values of cu are 
essentially unchanged from their x-^ = value. 

Next, we compare the thermally averaged scattering rate to the single-nucleon 
scattering rate F^g^^, i.e., the scattering rate expected if there were no multiple- 
scattering suppression. Using ttItt = 0, which maximizes S^^, and approximating 
1 + f^- 1 for thermal pions we have St^{u;) — {AaT^/m\) XnC'^oJ where = 
Cl,o + CXi + 0. Then, 

2 /oAATv/r„Trx 2 



Tscat _ 10^.7-^ ^0 17 ( T X / SOOMeVX ^ 



scat '""'AT 

where we have set XnCl/{{Sv + 35^) /4) = 1. By using the single nucleon result, we 
have probably overestimated the contribution to Fgcat from nucleons; however, any 
multiple scattering suppression of Snn should also affect S^. Further, we have over- 
estimated the contibution from pions by taking m-j^ = 0. Therefore, the interaction 
of thermal pions with nucleons yields only a modest correction to the scattering rate 
for the conditions pertaining to a SN core. 

Now turn to the possibility of a pion condensate which implies that /lyr — 
IJ"n — IJ'p — (^ttO where cu^q is the lowest energy value for the 7r~. The interaction with 
nucleons typically causes a dispersion relation for pions where this minimum occurs 
for a non- zero momentum kyro so that the 4-momentum describing the condensate is 
^ttO = {^-kO, kjro)- The occupation numbers for the condensate are given by fni^n) = 
(2%)^ Njj^S^ (k^^ — kyro), reducing the 7r~ phase space integration in S^'^ to N-,^ / {2uJt^o) 
and leaving a factor S{u!t^o — cu). 
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Most authors use a condensate with k^rO 7^ 0, based upon consideration of cold, 
degenerate nuclear matter characteristic of an older neutron star. It is conceivable, 
however, that thermal effects may modify the pion dispersion relation so that ^^^{k.;^) 
has its minimum at kT^o = 0. Wc will first consider such a zero-momentum condensate. 
The contraction Sl^^N^jj is identical to that for thermal pions except that we have to 
set both uIt^ = u!.,^q and m^r = i^nO- Then we find 

S^ico) = IGn^a^ (cXq + 3Cf^^o) fc^Ko + ^) + ^n(^Ko - ^)] ■ (5.14) 

We have used the fact that 1 + ^ /tt for the condensate. Then Xp/Xn = e~'^'^°/-^, 
and detailed balance is satisfied. Curiously, the axial and vector charges have ex- 
changed roles from the single nucleon fact which is related to the 75 in the 
pion coupling to nucleons. 

It is easy to compare scattering and pair absorption for l.h. neutrinos by com- 
paring Fscat with Fpair in Eq. (4.7) at a; = uJt^q/T. The pion contribution to the pair 
process exceeds that to scattering for ui^^q > GAT. Equally, we may compare the 
energy-loss rate in r.h. neutrinos from Tr~p nui^uji with that from i'i7i~p — > nuji 
or 1/1,71 T:~pvji. Here, the comparison is between Fgcat with Fpair of Eq. (5.5). Pair 
emission exceeds spin-flip scattering for uj^^q ^ 5.5 T. Therefore, in both cases scatter- 
ing is more important than the pair processes because for T fti 50 MeV a condensate 
with uIt^q ^ 5 T is highly unlikely. 

Next, we compare the pion-induced scattering rate with the single-nucleon one. 
In analogy with Eq. (5.13) we now find 

J- scat ^N-^ 

with X = uj-j^q/T. For example, taking the number density of pions to be N^^ = O.INb, 
ruN = 800 MeV, T = 50 MeV, and c^^rO = 100 MeV we get E^^at/r^at ~ 0-2- Even 
for small values of cUttq it is difficult to get pion-induced scattering to dominate over 
that from single nucleons. Further, as uJt^q is decreased the result will become more 
sensitive to the scheme used to include multiple-scattering effects for soft processes. 
To conclude, a zero-momentum condensate is unlikely to dominate over single nucleon 
scattering for neutrino transport or the emission of r.h. neutrinos. 

Finally, we consider finite momentum condensates. As long as (wtto, kTro) is 
time-like, we may always go to a Lorentz frame with kv^o = at the expense of 
exact isotropy of the nucleon and neutrino distributions. If |k7ro|/<^7rO is not large, 
then in the new frame the fluid momentum will be small, number densities will be of 
order their values in the fluid rest frame, and we may still use the long- wavelength 
approximation. Apart from an overall factor of order unity involving cUj^q and kTro, 
there will be no dramatic changes from the case of a zero-momentum condensate. 

However, when the condensate dips below the light-cone so that {wnO, ^no) is 
space-like, cu^g < |k7r0p) a different approach is required. In this case the denominator 
[(p+^7r/2) -/ctt]^ in Eq. (5.11) can become zero, corresponding to on-shell intermediate 
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nucleoli states. Put differently, the process pTr~ ^ nis now possible without any other 
particles involved. However, because the pions have one fixed momentum, we may 
consider a Lorentz frame where lo'^^q — 0. There, the pion condensate looks like a static 
pion field with a fixed wave-number k^Q. The nuclcons, then, should be described as 
Bloch waves in this periodic potential, i.e., the nucleon quasi-particles will be certain 
superpositions of n and p, involving spatial Fourier components of typical nucleon 
momenta and of k^g. 

The NC scattering of neutrinos from these quasi-particles in the non-rclativistic 
and long- wavelength limits should be very similar to the scattering off quasi-free single 
nucleons. Therefore, in the non-degenerate limit it should be given essentially by the 
number density of baryons times a typical weak cross section on a nucleon. We do 
not expect an anomalously enhanced scattering rate. (We note that in old neutron 
stars a pion condensate leads to a strong enhancement because the pion momentum 
is available to conserve momentum and thus, processes can go which otherwise are 
suppressed by degeneracy effects. In our case, there are no barriers from momentum 
conservation to overcome, and scattering can occur with full strength anyway.) 

To summarize, we arrive at two conclusions. First, pions would always enhance 
the scattering rates more than they would enhance pair processes, whether or not 
they are in a condensate. For r.h. neutrino emission that was considered in Ref. 8 
this means that its effect would be most important for the spin-flip scattering channel 
nvji. Second, however, we find that typically the effect of pions will not 
yield an anomalous enhancement of the scattering rate relative to the single nucleon 
case. The overall uncertainty of the neutrino interaction rates at high densities ap- 
pears to be much larger than the uncertainty introduced by the question of whether 
or not there is a pion condensate. 

We must comment that our conclusions differ somewhat from Mayle, et al. [10], 
who find that a pion condensate can have a dramatic effect, especially for the emission 
of r.h. neutrinos. There are two essential differences. First, they estimate the emission 
of r.h. neutrinos by multiplying the emissivity per thermal pion times the number of 
pions in the condensate, which we believe significantly overestimates the emissivity. 
For example, we can compare the emissivity, Q, per pion by using Eqs. (5.12) or 
(5.14) for S'tt, for thermal pions or for a kv^o = condensate, respectively. Plugging 
into Eq. (5.4) and using Fpair, the ratio of the emissivity per thermal tt divided by 
the emissivity per condensate tt is. 



Unless the condensate is "hard" , use of the thermal pion formula overestimates the 
emissivity by an order of magnitude; but, if lOt^q > T then Eq. (5.15) shows that 
the rates will be small compared to those from single nucleons. Further, as argued 
above, for a kTro ^ condensate we expect little enhancement over the single nucleon 
emissivity, if one uses an appropriate set of Bloch states to describe the nucleons. 

We have conducted our entire discussion of the pion-induced effects as if we 



(Q/iV^)th 

((5/A^7r)cond 



30C2 



(5.16) 



(Cl,o + 3C^,o)^§(l + e-^°)' 
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were in a dilute medium where different contributions to tlie emissivities add linearly. 
The pion interactions, however, are just another contribution to the nucleon spin 
fluctuations and so, they essentially add to the width of Sa{co) without adding to the 
overall strength, at least as long as we may treat the normalization in Eq. (3.3) as an 
upper bound. If we take our conclusions concerning S'tvtv and S'j^j,^ as a guide, it is 
quite possible that a large number of pions could actually have the effect of decreasing 
the effective emissivities of axions or r.h. neutrinos. 

The second distinction is that Mayle, et al. [10] stress that the presence of 
charge in a pion condensate will reduce the electron degeneracy, which in turn results 
in the release of entropy and heating of the core material. The resulting increase in 
temperature increases all emission rates. In this picture, the presence of pions affects 
the emission rates indirectly through their impact on the equation of state. This 
aspect of Maylc et a/.'s work is precisely in the spirit of what we advocate: When 
adding novel physics to the model one must be careful to include all ramifications, 
not just focusing on a single aspect of the new phenomena. In the present context, 
if one adds a pion condensate to the model to increase the emission of r.h. neutrinos 
then one must also allow for a change in the transport properties of the ordinary 
l.h. neutrinos. 



VI. DISCUSSION AND SUMMARY 

We have studied weak NC processes in the environment of the newly born neutron 
star in a SN core. They are crucial for understanding the cooling of the core for the 
first tens of seconds after collapse, the time during which the neutrino fiux from a 
galactic SN would be observable in a large underground detector. The most important 
application is the role NCs play in ordinary neutrino transport of lepton number and 
heat. In addition, they also couple weakly to a variety of hypothetical particles, 
which would provide for an anomalous loss of heat from the core by particle emission. 
Constraints on exotic particle properties can be derived from a comparison between 
their calculated effect on the neutrino fiux from the "surface" of the star with the 
SN 1987A observations. 

Understanding these transport and emission processes is made difficult because 
of the breakdown of perturbation theory, the usual tool for studying particle interac- 
tions. Nonetheless, it is possible to relate different aspects of a particular problem, 
e.g. the transport of lepton number or energy, to a common medium structure func- 
tion, S{u>), where uj is the energy transfer to the medium. The rate for any weak 
process can then be put into the form / du F{lu)S{u), where F is specific to that 
process. The difficulty with perturbation theory is confined to the function S{u!), so 
that once a particular S{u!) has been determined, the rates for all related processes 
can be calculated in a consistent fashion. Of course, different processes have F's 
which weigh the frequency dependence of S in different ways, so there remains some 
sensitivity to the shape of S{u!). Within this framework we have compared neutrino 
scattering and pair absorption and emission as they apply to the transport of ordinary 
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neutrinos. We have also studied how the same functions S{u!) relate to the emission 
of r.h. neutrinos and axions. 

In the limit that the medium is dominated by non-relativistic nucleons, the 
medium response consists of the sum of vector and axial vector response functions, 
Sv{^^) and Sa{uj) respectively. In the context of perturbative calculations of S, we 
have discussed two contributions: one from quasi-free nucleons So6{u!), and one from 
interactions in the presence of two nucleons Snn{u). In the non-relativistic long- 
wavelength limit the former is a delta function, whereas the latter naively has a soft 
divergence Snn{uj) ~ uj~'^. In § III we argued that, based on physical arguments and 
explicit calculation, Snn{u!) contributes only to Sa, i-e. the vector response is well 
approximated by quasi free nucleons. 

The axial response on the other hand requires a full calculation of the spin 
fluctuations in the medium, and is not well determined. Nucleon-nucleon interactions 
randomize the nuclear spin at a rate comparable to the collision frequency, which has 
a two- fold effect. They soften SaS{w) to something like a Lorentzian with a width of 
about the nucleon-nucleon interaction rate Tint and, by the same token, they regulate 
the divergence of Snn{u!) at = 0. It is possible to choose the width of the Lorentzian 
so that it reproduces the quasi-free nucleon 5-function result in the limit of a very 
dilute medium, while at the same time for large values of uu the bremsstrahlung result 
for 5'^7v(cc;) is obtained. Although this formulation is appealing (it replaces the one 
and two nucleon results by a single function), if extended into the regime of large 
collision rates it implies a very strong suppression to the axial response function and, 
as discussed in § III, the internal consistency of the approach may be questioned, as 
well as its justification based on a comparison with laboratory data. 

We conclude that there is no unique way to extrapolate the Snn to high den- 
sities. One can obtain anything from an enhanced to a much suppressed neutrino 
scattering rate. We stress that this inadequacy of perturbation theory cannot be 
ignored. We are interested in the evolution of a neutron star for, say, the first ten sec- 
onds after the collapse of the progenitor star. During this phase, all numerical models 
show maximum temperatures around T — 50 MeV. Even though large parts of the 
core may be cooler, the neutrino opacity increases with temperature, and so the hot 
regions should act as a bottleneck to neutrino transport. There, the nuclear medium 
is on the verge of degeneracy, but still essentially non-degenerate. If the medium were 
sufficiently dilute, then even for strong interactions one could proceed perturbatively 
with some confidence, while if the degeneracy were much higher only nucleons near 
their Fermi surfaces participate, in effect rendering the medium more dihitc. In the 
middle, perturbation theory fares worst. Still, because one can calculate reasonable 
answers both for very high and for very low degeneracies, we find it difficult to believe 
that the neutrino scattering rate should be very different from its naive lowest-order 
result in the intermediate range. In particular, an anomalously enhanced scattering 
rate seems rather unlikely, especially as it would require a violation of the upper limit 
in Eq. (3.3), which is based on the presumption that the different momentum, spin, 
and isospin parts of the nucleonic wave function are uncorrelated. 

We speculate that the overall shape of S{u;) given by the low-density approxi- 
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mation will not be radically modified at high densities. In this case scattering rates 
definitely exceed the rates for pair processes. This conclusion applies to the transport 
of l.h. neutrinos as well as to the emission of r.h. ones. Moreover, extreme modifica- 
tions to the shape of Sa seem likely to suppress the total interaction rate for both 
scattering and pair processes. Since the vector response is essentially unmodified and 
contributes only to scattering, it seems unlikely that any reasonable form for Sa will 
make pair processes dominate over scattering. This conclusion holds equally for the 
scattering of l.h. neutrinos and the emission of r.h. neutrinos. 

We caution, however, that there seems to be a trend that as the medium be- 
comes less dilute the response function gets suppressed for low cu while it increases on 
the wings. We have assumed that for cu above a typical nucleon-nucleon interaction 
rate, S is still reasonably represented by a perturbative result. If this were not the 
case, the true shape of S{u) could be so "hard" that pair processes could be important 
for neutrino transport after all. Equally, if the medium had strongly excited collective 
modes, their decay into neutrino pairs could be stronger than their contribution to 
scattering, depending on the nature of their dispersion relation. In this light we note 
that Iwamoto and Pethick [6] have discussed the importance of sound waves with a 
strong NC coupling, while Haensel and Jerzak [25] have shown that in degenerate 
matter quasi-free nucleons still seem to make the dominant contribution. Further, 
in the explicit case of a time-like 7r~ condensate we find that pair emission does not 
dominate. 

Given the uncertainty in the basic transport rates, how is one to understand the 
good agreement between the SN 1987A neutrino signature with theoretical expecta- 
tions? One answer is that other parameters may have been adjusted to achieve good 
agreement with the SN 1987A neutrino signal, hiding the sensitivity to the neutrino 
transport. Another possibility is that somehow the cooling curve is rather insensi- 
tive to the details of neutrino transport. Or, perhaps the transport in the numerical 
models, just by chance, is close to the real thing. For example, if the axial response 
were eliminated entirely, the vector resoonse would still contribute 1/5 of the naive 
single nucleon scattering rate. Perhaps even a fifth of the standard scattering rates 
is sufficient to reproduce the observed cooling timescale of SN 1987A. In this light, 
we feel it is time to undertake a systematic survey of how variations in the neutrino 
diffusion affect the cooling of the core on timescales from a second to a minute. This 
certainly seems important in light of the efforts to prepare experiments in advance of 
the next galactic supernova neutrino burst. 

We have also commented on the derivation of bounds on exotic particle interac- 
tions. In such calculations, the non-standard processes should be treated on the same 
footing with the ordinary neutrino transport, i.e., all NC type processes should be 
based on the same function S. Any regulation applied to keep the ordinary scatter- 
ing rates bounded must be applied consistently to exotic emission rates also, which 
essentially has the effect of decreasing them. Further, new processes proposed to in- 
crease emission rates must be studied for their effects on neutrino transport as well. 
In the end, only calculations where the transport and emission are performed with 
the same approximations for S{u!) may be used for quantitatively precise constraints. 
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and even then only after allowing for different choices of S{u!) and variation of all 
other unknown parameters. 

Finally, we have evaluated the contribution to the response function from the 

interaction of single pions with nuclcons, Sty{uj). Contrary to previous speculations 
in the literature, we find no indication for an anomalously large contribution by this 
latter process for the conditions relevant for the cooling of a SN core. 
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APPENDIX A: EFFECTIVE NUCLEON WEAK CHARGES 



If the nucleon weak current is written in the form Eq. (2.3) the vector charge is given 
by Cv = T3 — Asm^ dwQ) where = +1 for protons and —1 for neutrons and Q is 
the nucleon electric charge in units of e. Therefore, 

C{> = -1 and C^ = 0.07, (A.l) 

where we have used sin^ 9w = 0.2325 for the weak mixing angle. 

There remain problems, however, with the value of the axial charge. If axial cur- 
rents were conserved we would have Ca — T3, but axial charge conservation is known 
to be violated by the strong interactions, leading to a charged-current axial coupling 
of bare nucleons of Ca,cc.Q = 1-26. In large nuclei this value is suppressed somewhat, 
and the commonly used value [26] for nuclear matter calculations is Ca,cc = 1-0. 
Based on the naive quark model it was common practice until recently to take the 
NC axial charges to be isovector in character and so, Ca,nc = 1-0 rs was a reasonable 
choice. It is now realized, however, that the neutral axial current has an isoscalar 
piece as well, which is associated with the polarization of the strange quark sea of 
nucleons. 

The axial charge of nucleon i can be written as a sum over the contributions 
from different quarks, = Am' — Ad' — As*. The contributions from u and d 
quarks include both valence and sea contributions, whereas the s quark only has a 
sea contribution. It is expected that the heavy quarks do not contribute significantly. 
Isospin symmetry dictates that Au^ — AdP and Ad'^ — AvP. Charged current 
processes determine Am^ — Ad^ = 1.26 for bare nucleons. Scattering experiments of 
polarized muons on polarized hadronic targets by the EMC [27] and SMC [28] groups 
at CERN, and by polarized electrons at SLAC (E142 [29]) give a range of values for 
AsP between and —0.2. In a review of the phenomenolgy and theory, EUis and 
Karliner [30] suggest As = -0.11 ± 0.04, leading to C^ q = 1.37 and C\q = -1.15 
in vacuum. For nuclear matter we find 

C3 = -0.91 and = 1.09 (A.2) 

if Au, Ad, and As are suppressed by an equal factor 1/1.26. We do not give un- 
certainties because of the many possible sources of systematic errors, both in the 
interpretation of the laboratory results and in the extrapolation to nuclear matter. 

For the scattering of neutrinos on quasi-free nucleons the relevant combination 
of coupling constants is = j{Cy + 3C^) for which we find 

_ ( 0.89 for protons (1.41 in vacuum), ^. 
\ 0.87 for neutrons (1.24 in vaccum). ^ ' ' 

Therefore, the scattering rate is suppressed by about a factor of 0.63 for protons and 
0.70 for neutrons relative to their vacuum values. 
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FIGURE CAPTIONS 



FIGURE 1 

The functions Fgcat and Fpair defined in Eqs. (4.5 and 4.6), and the analytic approx- 
imations defined in Eq. (4.7). The results are normalized to Fq = GpN^T'^- 

FIGURE 2 

Comparison of scattering and pair absorption rates, Eq. (4.4), calculated using the 
F's in Eq. (4.7), for different choices of S{u;). Horizontal hnes are single- nucleon 
scattering rates. Dotted curves use Snn from Eq. (3.5) with a lower cutoff ci^nim while 
the dashed curves are for from Eq. (3.8), also cut off below cjmin- The parameters 
used are p = po = 3 X 10^^ gm cm'^, T = 50MeV, Xn = 0.7, = 800 MeV, and 
the vector and axial charges given in Eqs. (A.l) and (A. 2), respectively. 

FIGURE 3 

Comparison of scattering and pair emission contributions to the emission of r.h. neu- 
trinos from Eqs. (5.4) and (5.5), for the same choices of S{u;) and parameters as in 
Fig. 2. 
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